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Abstract 

We illustrate how to compute local risk minimization (LRM) of call op¬ 
tions for exponential Levy models. We have previously obtained a repre¬ 
sentation of LRM for call options; here we transform it into a form that al¬ 
lows use of fhe fast Fourier transform mefhod suggested by Carr & Madan. 

In particular, we consider Merton jump-diffusion models and variance gamma 
models as concrete applications. 
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1 Introduction 

Local risk minimization (LRM), which has more than 20 years' history, is a 
well-known hedging method for contingent claims in incomplete markets. Al¬ 
though its theoretical aspects have been very well studied, corresponding com¬ 
putational methods have yet to be thoroughly developed. This paper aims to il¬ 
lustrate how to numerically calculate LRM for call options in exponential Levy 
models. To our knowledge, this contribution is the first to address this sub¬ 
ject. In Arai & Suzuki |T], we obtained a representation of LRM for call options 
by using Malliavin calculus for Levy processes based on the canonical Levy 
space. Here we transform that result into a form that allows the fast Fourier 
transform method suggested by Carr & Madan to be applied. In particular, 
Merton jump-diffusion and variance gamma models, being common classes 
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of exponential Levy models, are discussed as concrete applications of our ap¬ 
proach. 

Consider a financial market composed of one risk-free asset and one risky 
asset with finite time horizon T > 0. For simplicity, we assume that the mar¬ 
ket's interest rate is zero, that is, the price of the risk-free asset is 1 at all times. 
The fluctuation of the risky asset is assumed to be described by an exponential 
Levy process S on a complete probability space (O, T, P)J3 described by 

Sf := So exp -F nWf-h j xN([0, t],dx)| forte [0,T], 

where Sq > 0, e M, u > 0, and Rq := R \ {0}. Here W is a one-dimensional 
Brownian motion and N is the compensated version of a Poisson random mea¬ 
sure N. Denoting the Levy measure of N by v, we have N ([0, t], A) = N ([0, t], A) 
tv {A) for any t e [0, T] and A e S(Ro). Moreover, S is also a solution to the 
stochastic differential equation 


dSt = St— 


}i^dt + t7dWt+ [ {e^ — l)N{dt,dx) , 
JtRo 


where •.= }i + \o''^ + — 1 — x)v{dx). Without loss of generality, we 

may assume that Sq = 1 for simplicity. Now, defining Lt := logSf for all 
fe [0,T], we obtain a Levy process L. Moreover, dMt := Sf- [udWt-|- — 

l)N{dt,dx)'\ is the martingale part of S. 

Our focus is the development of a computational method for LRM with re¬ 
spect to a call option {Sj — K)^ with strike price K > 0. We do not review the 
definition of LRM in this paper; for details, see Schweizer ( fT6l , IflTl l. We first 
briefly introduce the explicit LRM representation of such options in exponen¬ 
tial Levy models given in [IJ. 

Define the minimal martingale measure P* as an equivalent martingale 
measure under which any square-integrable P-martingale orthogonal to M re¬ 
mains a martingale. Its density is then given by 

^ = exp I Wt + log(l - 0x)N( [0, T], dx) + tJ^ e,v{dx) | , 

where 

++-'^Vv{dy) 

for X G Rq. In the development of our approach, we rely on the following: 

Assumption 1.1. 1. J^{\x\\/x^)v{dx) < oo, and — l)’^v(dx) < oofor 

n = 2,4. 

^ (fl, A,P) is taken as the product of a one-dimensional Wiener space and the canonical Levy 
space for N. Moreover, we take F = {At}t£[ o,rl the completed canonical filtration for P. For 
more details on the canonical Levy space, see ll9l and (T]. 
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2. 0>^^ > -cr^ - - l)^v{dx). 

The first condition ensures that and 9x are well defined, the square in- 

tegrability of L, and the finiteness of — l)^v(dx) for n = 1,3. The sec¬ 

ond guarantees that 9x < 1 for any x e Mg- Moreover, by the Girsanov theo¬ 
rem, Vlf* := Wf -|- and N^*{[0,t],dx) := 9xv{dx)t + N{[0,t],dx) are a P*- 
Brownian motion and the compensated Poisson random measure of N under 
P*, respectively. We can then rewrite Lf asLf = fi*t + o'Wf'* + xN^* {[0,t],dx), 
where }i* := -|- 1)(1 — 9x)v{dx). Note that L is a Levy pro¬ 

cess even under P*, with Levy measure given by v^*(dx) := (1 — 9x)v{dx). 
The LRM will be given as a predictable process LRMt, which represents the 
number of units of the risky asset the investor holds at time t. First, we define 


h ■= [I{Sr>K}ST I J^t -], (1) 

h := [ Ep* [{Sre^ - K)+- {St - K)+\ - l)v{dx). (2) 

JRo 

Our explicit representation of LRM for call option (Sp — is then as follows: 
Proposition 1.2 (Proposition 4.6 of U). For any K > 0 and t e [0, T], 


LRMt 


_+ h _ 

St- 


(3) 


Remark 1.3. 1. The assumption {e^ — l)'^v{dx) < ca is imposed in Proposi¬ 

tion 4.6 o/ifit. 

2. If the interest rate of our market is instead r > 0, then Q becomes 


LRMt = 


_ cr^h + h _ 

+ -l)^v{dx))' 


and P* is rewritten with L and 9x becoming ?, r, „ X-, , ■ ^ 

S X 6 ir2+J^{ey-l)2v{dy) 

and , ^ 1 - —ri7^vr4j-^/ respectively. Moreover, the second condition in As- 

sumption \1.1\ would be revised to 0 > — r > — l)^v(dx). 

That is, a nonzero r requires only that we replace y with y — r and multiply the 
the expression for LRMt by which means that we can easily generalize 

results for the r = 0 case to those for r > 0. For simplicity, in this paper we 
treat only the case r = 0. 

From the point of view of Proposition [L2l we have to calculate conditional 
expectations of functionals of Sj under P* in order to calculate LRMt numeri¬ 
cally. However, there does not appear to be any straightforward way to specify 
the probability density function of Sj (or equivalently Lp) under P*. Instead, 
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since L is a Levy process, it may be comparatively easy to specify its charac¬ 
teristic function under P*. Hence, a numerical method based on the Fourier 
transform is appropriafe for computing LRM. Moreover, Carr & Madan [4J in¬ 
troduced a numerical method for valuing options based on fhe fasf Fourier 
transform (FFT). We take advantage of fhis fo develop a numerical method for 
LRM. To fhis end, we induce integral expressions for Ii and I 2 in terms of the 
characteristic function of Lr-t under P* and recasf them into a form thaf al¬ 
lows fhe Carr-Madan approach to be applied. In particular, I 2 will be given as 
a linear combination of Fourier fransforms. 

In this paper, we consider two concrete exponential Levy processes for L. 
The first is a jump-diffusion process as introduced by Merton ml This con¬ 
sists of a Brownian mofion and compound Poisson jumps with normally dis¬ 
tributed jump sizes. The second is a variance gamma process, which is a Levy 
process with infinitely many jumps in any finite time interval and no Brownian 
component. This was introduced by IIT2I and can be defined as a time-changed 
Brownian motion. Many papers (e.g., [4j, IITSl l have studied it in the context of 
asset prices. Schoutens IITSi provides more details on these two Levy processes 
and more examples of exponential Levy models. 

There is great deal of literature on numerical experiments related to LRM 
(e.g., f3j, fZj, fSj, IITOl . lITTl . | 21] ), but to our knowledge, ours is the first at¬ 
tempt to develop an FFT-based numerical LRM scheme for exponential Levy 
models. Kelani & Quittard-Pinon [91 studied an optimal hedging strategy that 
they call 6-hedging, which is similar to but different from LRM, for exponential 
Levy models, and adopted a Fourier transform approach separate from Carr 
& Madan U's method. As an important difference, they assumed that S is a 
martingale under the underlying probability measure. In contrast, we do not 
make this assumption. We therefore need to treat S under P*, that is, calculate 
conditional expectations of functionals of S under P*. However, the structure 
of S is no longer preserved under a change of measure. For example, when L 
is a variance gamma process under P, it is not so under P*. Thus, our setting 
is more challenging but also more natural. 

The rest of fhis paper is organized as follows: An introductory review of 
the Carr-Madan approach is given in Subsection 2.1, and the integral repre¬ 
sentations of Ii and I 2 are presenfed in Subsection 2.2. Merton jump-diffusion 
models are examined in Section 3, which starts with mathematical preliminar¬ 
ies and proceeds to numerical results. Section 4 is similarly devoted to variance 
gamma models. 

^ Merton m also suggested a hedging method for these models, but this is is different from 
LRM. For additional details, see Section 10.1 of (5|. 
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2 Preliminaries 


2.1 Numerical method 


We briefly review the Carr-Madan approach, which is an FFT-based numerical 
approach for option pricing. The FFT, introduced by |[6l, is a numerical method 
for computing a discrete Fourier transform given by 

N-l 

F{1) := (4) 

j=0 


for / = 0,..., N — 1, where {^y}y=o,...,N-i is a sequence on IR and where N is 
typically a power of 2. The FFT requires only 0(Nlog2N) arithmetic opera¬ 
tions, as compared with the usual Fourier transform method's O(N^). 

The aim of the Carr-Madan approach is efficient calculation of E[(S 7 ’ — 
X)'*’] when S is a P-martingale. Recall that we are considering only the case in 
which the interest rate is zero. Denoting k := logK and C(fc) := E[(S 7 ’ — 
we have 


C{k) = - r dv 

n Jo liv — ia)\t[v — toi) + 1\ 


(5) 


for a > 0 with E[Sj“'’^] < oo, where cp is the characteristic function of Lj. 
Note that the right-hand side of ll5]l is independent of the choice of a. Now, we 
denote tp{z) := z G C. Using the trapezoidal rule, we can therefore 

approximate C{k) as 




( 6 ) 


j=0 


where N represents the number of grid points and rj > 0 is the distance be¬ 
tween adjacent grid points. The right-hand side of ([bll corresponds to the inte¬ 
gral in ||5j over the interval [0, Nr/], so we need to specify N and t] such that 


1 POO 

- / - ioiidv 

71 JNij 


< e 


(7) 


for a sufficiently small value e > 0, which represents the allowable error. By 
incorporating Simpson's rule weightings, we may rewrite ijb) as 


C{k) « i "y: - k)|(3 + (-ly+i - 5j), 


M 

where <5, is the Kronecker delta function. We define 


o—cik N—1 


j=o ^ 


71 
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for I = 0,...,N — 1, which is a discrete Fourier transform as given in l(4]l. This 
yields 


C{k) 


ZTT 


So long as we take rj so that |A:| < n/rj, we can employ the FFT to compute 
C{k). 


2.2 Integral representations 

We next induce integral expressions for fj and I 2 , defined in |[T]| and (|2]l, and 
evolve them so that the Carr-Madan approach is available. Recall that As¬ 
sumption [id] applies throughout. As can be seen from Subsection 2.1, if 7i and 
I 2 are represented in the same form as l|5| we can compute them by means of the 
Carr-Madan approach. Because the conditional expectations appearing in h 
and I 2 are under P*, the functions corresponding to in (|5]l should include the 
characteristic function of Lj_t under P*, denoted by (pj_t{z) := Ep* 
for z e C. 

First, we induce an integral representation for Jj (=Ep« [l{Sj.>jc}Sr | ^t-]) 
with (pT-t by using Proposition 2 from Il20l : 

Proposition 2d. For K > 0, 

1 ®+i 

^tp*['^{St>k} ■ St \ :Ft-] = - ( 8 ) 

for all t G [0, T] and a. G (1,2]. Note that the right-hand side is independent of the 
choice of OL. 

Proof. Define G{x) := 1{3:>k} ■ x, g(x) := G(e^) for any x G R, and g{z) : = 
/r for any z G C. We employ one lemma: 

Lemma 2.2. Let L'be an independent copy of L. Then, LT_f + N- = Lj for all 
t G [0, T], where A B means that A = B in law for P*. 

Proof of Lemma 12.21 Proposition 1.7 of t2j implies that P*(Lf_ = Lf) = 1. 

F*-d 

Therefore, Lf = Lf_. Because Levy processes have independent and station- 

'P*-d 

ary increments, we have Lj = Lj — Lt Lf = L'j_f -|- Lf. □ 

Returning to the proof of Proposition 12.II from Lemma I2l2l we have 

Ep*[l{Sr>ic} ■ St I d^t-] = Ep* [G(S7’) I J'f-j = Ep»[g(Lp_, -h Lf_) | Ft-] 

= / g{x + Lt-)p{dx), 

^]R 

where p(A) := P*(Lj_j G A) for any A G S(R). By (22)-(25) in the proof of 
Proposition 2 of 1201 , if any a G (1,2] satisfies the conditions that 
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(a) g{x)e has finite variation on IR, 

(b) g{x)e-^^ G Li(R), 

(c) Ep* < 00 , and 

then 

f g{x + Lt-)p{dx) = ^ f g{v + ia)cpT-t{-'v-ioc)SfZ'^dv 

for a G (1,2], which is independent of the choice of a. As a result, under 
conditions (a)-(d), we have 


Eit 


■ [l{Sr>ic} ■ Sr I ^ J^g{v + m)4)T-t{-v-m)S^_'^dv 

= — g{-v + ia.)^T-t{'v-ioi)Sf:!^'^dv 


1 r°° 

— - —^(pT-t{v-m)Sf+'^dv. 

n Jo x — \ + iv ^ 


We need only to confirm that conditions (a)-(d) hold. Conditions (a) and (b) 
are obvious. To demonstrate condition (c), it suffices to show Sp_t G L2(P*) 
for any f G [0, Tj. Note that we have 


I {e^-lfv^*(dx) 

JKn 


[ {e^-l)^v{dx) + ^ - j :— ^ , . [ {e^-lfv{dx) 


< 00 . 


Because S is a solution to dSt = St-{crdWf* + /^(e^ — 1)N^*{dt,dx)), Theo¬ 
rem 117 of fTS i implies that sup^g p j] | St | G (P*). 

Next, we show condition (d). Note that 


— ioc) 


= Ei; 


exp < {iv + ol) 


f xN^'’ {[0,7 - t],dx) 

./Ro 


(9) 


For the right-hand side, we have 


Ep* 


exp 


{iv + tx) f xN^*{[0,T — t],dx) 
JlRg 


^ Ep* 


exp •] a [ xN^*{[0,T — t],dx) 

Jro 


< 00 , 


( 10 ) 
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because 

Ep. 


— Ep* 


exp < a. 


^*{T-t)+ aWflf + I xN^* ([0, T - f], dx) 
JRn 




ilM-t],dx) 


Ep* = exp 

obtain 


Ep* [exp{ {iv + a)(7Wp_f }] 
As a result, we have from (l9ll-lfTTll 


I j 0 c^a'^{T — t) and Ep* < 00 . In addition, we 

(a^ — v'^)cr^(T — t) 


f m-i 

JlR 1 


+ P 


-dv < C 


It 


= exp 


exp 


(11) 


+ \v 


2 , , 

--- -V > dv < CO 


for some C > 0. This completes the proof of Froposition l2.ll 
We evolve l|8]l into the same form as lO as follows: 

1 f°° 


□ 


h = Ep*[l{sr>j<c} ■ St I J^t-] = - 


(pT-t{v — m)Sl^''^dv 


x — \ + iv‘ 

rOO 

/ ( 12 ) 
Jo 


d>T 

where k := logK and !/^i(z) := ^ /z-i for z e C. Thus, we can compute h 
with the FFT based on Subsection 2.1. 

We turn next to J 2 (= Ep* [(Sje^ — K)+— (Sj — ff)+ | — 'J)yidx)). 

First, we have the following integral representation: 

Proposition 2.3. For any K > 0, 


Ep* [(St-X)+ 


= 


1 r 

n Jo 


K 


-iv—oc-^l 


cpT-t{v-m)S^+^^ 
(a — 1 + iv){oi + iv) 


dv 


(13) 


for any t e [0, T] and any a e (1,2]. Note that the right-hand side is independent of 
the choice of x. 

Proof. We can see this in the same manner as Proposition l2.1l but with G{x) = 
{x-K)+. □ 

Note that dlSll coincides with ||5]l, where a — 1 in dlSt corresponds to a in ||5j. 
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Denoting i/’ 2 (z) 


<Pt-M)S\L 

{iz—l)iz 


for z e C and ^ := v — ia, we have 


Ep. [(St-X)+ I Tt-] 


-I 

TL Jo 

TC Jo 


K 


K- 


—iv—a+l 


4>T-t{v - 

(a — 1 + iv){a. + iv) 

-i^+i 4>T-t{0sf_ 

POO 

K-^^+^^P2{0dv=:f{K). 


dv 


(14) 


Note that /(X) is computed with the FFT. Moreover, Fubini's theorem implies 


h 


f Ep. [(STe^-X)+- (St- X)+ | J'f_](e^ - l)v(dx) 

Jro 

f {e^f{e-’‘K)-f{K)}{e^-l)v{dx) 

JKq 

r (1 I 

/ {Ke~'‘)~'^^'^ip2iC)dv -/ K~'^+^ip2{0dv\ 

J IRq ^ 0 TT^O J 

I - l)X-'f+V2(Orf4 (e" - -^Mdx) 

J IRq 7T 7 0 J 

1 POO P 

- -l)v{dx)xp 2 { 0 dv, 

TC JQ j IRq 


l)v(dx) 


(15) 


which is the same form as 10, because the integrand of ifTSll is a function of 
However, we carmot compute (flSll numerically as it stands, because it is 
not possible to compute the integral — 1)(£^ — l)v(iix) directly. Thus, 

we need to make further model-dependent calculations. In Sections 3 and 4, 
respectively, we evolve dlSll into a linear combination of Fourier transforms for 
Merton jump-diffusion models and variance gamma models. 


Remark 2.4. Regarding LRMt, h, and I 2 as functions of St- and K, we have 
Ii{St-,K)/St- = Ii{l,K/St-)fori = l,2by®and ^,and 


LRMtiSt-,K) 


a^h{St-,K) + l2{St-,K) 
St_(cr2 + 4^(e^-l)2i/(dx)) 


u2fi(l,X/Sf-) + f2(l,X/St-) 


by 0. As a result, LRMt is given as a function of K/St- =: nit-, where nit- is 
called moneyness. Thus, we denote LRMt by LRMt{mt-). As a by-product of this, 
we can analyze jump impacts on LRM. If the process L has a jump with size y e IRq 
at time t, then the moneyness mt- changes into mt-e~y at the moment when the 
jump occurs. Thus, LRM also changes from LRMtjmt-) to LRMt{mt-e~y). We can 
regard the difference LRMt{mt-e~V) — LRMt(mt-) as a jump impact. In particular, 
LRMt{e~y) — LRMt(l) represents a jump impact when the option is at the money. 

Remark 2.5. Hereafter, we fix a E (1,2] arbitrarily. Moreover, we denote f := v — ia 
for i; e R, so we may regard ^ as a function ofv. 
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3 Merton Jump-Diffusion Models 

We consider the case in which L is given as a Merton jump-diffusion process, 
which consists of a diffusion component with volatility <7 > 0 and compound 
Poisson jumps with three parameters, m € R, <5 > 0, and 7 > 0. Note that 
7 represents the jump intensity and that the sizes of the jumps are distributed 
normally with mean m and variance Thus, its Levy measure v is given by 

= vis 

When it desirable to emphasize the parameters, we write v as ^[7, m, <i]. Note 
that the first condition of Assumption ll.ll is satisfied for any m € R, <5 > 0, and 
7 > 0 . In addition, the second condition is equivalent to 

0 > -H Y + 7 jexp (m + y'j - l-m\ 


and 


3(^2 


}i -\—-—h 7 s exp( 2 OT -h 23'^) — exp —ot>> 0 . 


<52 


We consider only the case in which the parameters satisfy Assumption ll.il 


3.1 Mathematical preliminaries 

Our aim here is threefold: (1) to give an analytic form for (pj_f(z) (:=Ept 

(2) to evolve ifTSt into a linear combination of three Fourier transforms; and (3) 

to give sufficient conditions for Ni] under which ((Tj holds for a given e > 0 . 

First, we provide an analytic form of (pr-t- To this end, we begin by calcu¬ 
lating v'*’*. 

Proposition 3.1. We have 


(dx) = v[{lh)j,m,S^]{dx)V —hyexp 

Proof. By Assumption [T tJ 0 > h > —1. Hence, 


f 2m -|- 3^ 


,m + 3^, 3^ 


{dx), 

( 1 ) 


v^* (dx) = (1 — 9x)v{dx) = (1 — h{e^ — l))v{dx) = (1 -|- h)v{dx) — he^v{dx). 
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Moreover, 


e^v{dx) = exp\x- \ dx 


y/lnS 

7 


exp 


2^2 


[x —(m + <52)]2 2m+ 5^ 


= V 


7 exp 


2m + <52 


2^2 


,m + S^,5^ {dx), 


dx 


from which l[D follows. 

Next, we calculate cJj-tiO fo’" ^ ^ [0/ '^]- 

Proposition 3.2. For any f £ [0, T] and u £ IR, with ^ := v — ia, 

u 2^2 


□ 


J>T-t{0 = exp|(r-f) 

= exp \{T-t) 




+ (1 + - 1 - imO 


—hje ^ S- — 1 — i^{m + 6^)\ 




Proof. We only have to show the first equality: 


‘pT-tiO — Ej 


exp } il, 


y* {T-t)+ crWflf + [ ([0, T - t\,dx) 

^Ro 


= exp {(T - t)iCy*}-Er. 


exp|!^ I xN^ ([0, T —f],dx)| 


= exp|(r-f) il,y*-^- + j^{e'^^-\-i^x)v^\dx) |. 

□ 

Second, we evolve lUSll . We define !^(z) := tp 2 (z) exp | —for z £ C 

and / {K) := ^ K~’^~^^{p{l,)dv. Remark that / is computed with the FFT as 

well as / defined in dlill . The following proposition demonstrates (ITSll . namely, 
I 2 is given by a linear combination of three Fourier transforms. 

Proposition 3.3. We have 

I Ep. [{Sje^ - K)+- {St - K)+\ Ft-] {e^ - l)v{dx) 

JRn 




^,-j{Ke-m-s^) _ ye’^f(^Ke-’^) + 7(1 - e’^+^ )f{K) ( 2 ) 


for any t £ [0, T]. 
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Proof. We calculate 

[ (e'f^-l)(e^-l)v(dx) 

./Ro 

= [ (e('f+i)^ - e'U + 1 _ e’‘)v{dx) 

./Rn 


/Ro 

= 7exp|(!C^ + l)m+ y(/^ + l)2| -7exp|f^m- l. 


,/7f + ^ > 


Hence, we obtain 


(flsll = — 

n Jo 

r\i p 00 c2 o x2 

ip2{0dv + y{l - e”^+^)f{K) 


= ye 


r/' 


fiKe-'”-^ ) - ye^f{Ke-^) + y{l - e^+^)f{K). 


□ 


Third, we provide sufficient conditions for the product Nt] under which ((7|l 
holds for a given allowable error e > 0. First of all, we determine an upper 
estimate for (pj-t- 

Proposition 3.4. We have 

\(pT-t{v - i(x)\ < Cl exp 
for any u e R, where 

Cl = exp|(T—t) ay* + ^^ + j {e'^^ — 1 — ax)v^ (dx) | 


= exp^T-t) 


cr^a^ 


ay* H— -h {l + h)y{e 




— 1 — am) 


2m+S^ 

—hye 2 


e(’”+^>+^ -l-a(m + (52) 


Proof. Proposition 13.21 implies that 

<pT-t{v-m) 

= exp|(T-h 

= exp((T-h 


a^{v — ictf 


i(v - m)iT - ^ - 1 - Uv - m)x)v^’ (dx)] \ 

2 ./Ro J J 

{iv + a)fi*- ^ r ^ ^ ^ (^dx)] \ 

2 /Ro J J 


= exp|(T—p/f 
X exp < (T — t) 


}i* + tPct^ f xv^'{dx)\\expl(T-t) f 
Jirq j j [ ./Ro 

(P(v^-a^) ^ r ^ 

JlR-n 


afi 


■+ / ( — l — ocx)v'^ (dx) 
Jrq 
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Noting that I exp {(T—1) <exp{(T—f) (dx)j, 

we have 


\(pT-tiv - ioi)\ < exp|(T-t) 


CCJi 


a^{v^ — cP-) 


+ [ (e“-l-ax)i/"’’(dx) I . 
JRq . J 


□ 


Propositions l3.5l and l3.6l below give sufficient conditions for Nr/ under which 
fi and U satisfy llZll for a given allowable error e > 0, respectively. 


Proposition 3.5. Let e > 0 and f e [0, T). When a > 0 satisjies 

1 


K r K 




we have 


< a, 


< e. 


(3) 


1 rco 1^—w—a+l 

- / -T— 

n Ja X - 1 -h w 

Proof. Noting that e~^ < x~^ for any x > 0, we have, by Proposition |331 

/•oo ]^—iv—a.+l 

/ - T—^(pT-t{v-m)S‘f+''’dv 

Ja 

1 /■“ 

n Ja |a — 1 + iv\ 

f K f°° 1 ^ crVfr f. 

\St- J Ja \x-\ + lv\ 


< 


K f K 

K / K 
n \St- 
K f K 

t- 

< e. 


Clri{^{T-t)} dv 

J a 


-At) -A. 


a V { 2 

-» .no 4^-5 


, K f 

rdV = — 


cA{T—tA ^\St-J (A{T—tAa 


Cl 


□ 


Proposition 3.6. Let e > 0 and f e [0, T). > 0 satisfies 

4Ci7X f K 
5ncA{T — tYe 

then 


|g(a+l)m+(.^+a+^),52 + |1 - 11 < 


1 pOQ p 

-/ (e'^^-l)(e^-l)i/(dx)!/^2(0^^^^ 

71 J a JIRq 


< £. 


(4) 

(5) 
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Proof. First, we estimate \^p 2 {I,)\dv. Noting that 


1 


1 

(ic-m 


{iv + a — l){iv + a) 


Proposition 13.41 implies 


I I 

Ja Ja 


(pT t{v-ioi)St 


i[v—ici) 


< 


{K-m 

4CiSf 


n4(r- f)2 
Hence, Proposition l3.3l implies that 


rc 

dv < CiSf_ / 
Ja 

r , 4CiSf 

/ v-^dv= 

Ja 




dv 


5(H(T — t)2fl5 ' 


L.H.S. of ID 


+ 


n 

7(1 - e^+T) 


Ja 71 Ja 

r K-‘^+^xp2{0dv 

Ja 


< 1 r L2m + iS^ 

71 Ja ' I 


e 2 


dv 


TT Ja 

.yh_e'«+T| ,1 

+ -1 / K-^+^\tp2(0\dv 

71 Ja 

- 5Tia^{T-t)^a5 [St-J 1 ' '/ 

< £ . 


□ 


3.2 Numerical results 

As seen in the previous subsection, substituting fl^ and l|2| for Ji and I 2 re¬ 
spectively, we can compute LRMt given in (|3l with the FFT. Note that we 
need Proposition 13.21 in order to calculate ipi, ip 2 , and {p. In this subsection, 
we provide numerical results for a Merton jump-diffusion model with param¬ 
eters T = 1, ji = —0.7, a = 0.2, 7 = 1, m = 0, and ^ = 1. Note that 
is given by —0.03, which satisfies the second condition of Assumption [O] In 

























3 MERTON JUMP-DIFFUSION MODELS 


15 


particular, we consider the following two cases: First, fixing the strike price K 
to 1, we compute LRMt for times t = 0,0.05,..., 0.95. Second, f is fixed to 0.5 
and we instead vary K from 1 fo 8 af sfeps of 0.25 and compute LRMq ^. Note 
that we take Lt- = 1 whatever the value of t is taken. Moreover, we choose 
N = 2^^, rj = 0.025, nad a = 1.75 as parameters related to the FFT. We have 
then Nrj = 409.6. For any parameter set mentioned above, both (|3) and (|4j 
are satisfied for e = 10“^. Figure [T] shows the results for these two cases. The 
computation time to obtain Fig. [ijb) was 0.59 s. Note that all numerical exper¬ 
iments in this paper were carried out using MATLAB (8.1.0.604 R2013a) on an 
Intel Core i7 3.4 GHz CPU with 16 GB 1333 MHz DDR3 memory. 
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(a) LRMt of a call option with strike price K = 1 and maturity T = 1 vs. time for a Merton 
jump-diffusion model with parameters = —0.7, a = 0.2, 7 = 1, m = 0, and S = 1. 


o 


8 


(b) LRMo ,5 for fho same Merton jump-diffusion model as (a) vs. strike price K. The verti¬ 
cal axis represents the value of 

Figure 1: Merton jump-diffusion model 
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4 Variance Gamma Models 


We now consider the case in which L is given as a variance gamma process. 
Note that L does not have a diffusion componenf. This means thaf <7 = 0, 
that is, fi vanishes. A variance gamma process, which has three parameters 
K > 0, m € R, and <5 > 0, is defined as a time-changed Brownian motion with 
volatility S, drift m, and subordinator Gf, where Gf is a gamma process with 
parameters (1 /k, 1/k). In summary, L is represented as 

Lf = ffiGf -|- for t E [0, T], 


where B is a one-dimensional standard Brownian motion. Moreover, the Levy 
measure of L is given by 


where 



‘=-55 




23'^ m 




m 


Note that C, G, and M are positive. To emphasize the parameters, we write v 
with parameters k , m, and <5 as v{dx) = v[k, m, <i] (dx). Moreover, by regarding 
G, G, and M as parameters, we may express v as v{dx) = addi¬ 

tion, we assume M > 4 in this section, which ensures that the first condition of 
Assumption 11.11 holds, by the following lemma: 

Lemma 4.1. When M > 4, — l)”v(dx) < ooforn = 2,4. 

Proof. For n = 2,4, we have 


< 00 , 


POO POO 

(e^ - l)”v(dx) <cj^ 

Cie^-lfvidx) < [\”{e-l)”v(dx) < Cie-lV <oo, 

Jo Jo 

j (e^ — l)”v(dx) < J {—x)"v{dx)<C j (—x)'’“^dx < oo, 
/ [e^ — l)”v(dx) < f v{dx) < C f e~^’‘dx 

J—oo J—oo J\ 


< 00 , 


because n — M<0, 0<e=^ — 1< x(e — 1) whenever x e [0,1], 1 -|- x < for 
any x 6 R, and < 1 if x < 0. □ 

Remark 4.2. We can generalize this lemma to |e=*= — l|“i/(dx) < oofor any a E 
[1,M). 




4 VARIANCE GAMMA MODELS 


18 


Because }i = xv(dx), (|2) below implies that the second condition of As¬ 
sumption [T^ can rewritten as 

, /(M-1)(G + 1)\ / MG \ 

°^(,(M-2)(G + 2)j ^ - °®(,(M-l)(G + l)j' 

which is equivalent to—3 < G — M < —1. 


4.1 Mathematical preliminaries 

The approach to variance gamma models is similar to that in Subsection 3.1. 
We begin by calculating of . 

Proposition 4.3. v^* (dx) = V(i+ft)c,G,M(^^) + ^ = 

Proof. By the same argument as Proposition 13.11 (dx) = (1 -|- h)v{dx) — 
he^v(dx). We have Avc,G.M(rf^) = i^ac,g.m(^^) for any A > 0, and 

e''vc,GMidx) = e^C(l{;,<o}e‘^^ -h 

= C(l{:,<o}e(^+^)^ -h = Vc,G+PM-l(rf^) 

because M — 1 > 0. □ 


Remark 4.4. For any A > 0, Av [k, m, <5] {dx) is a Levy measure corresponding to the 
variance gamma process with parameters k/ A, Am, and 3^/A- However, vc^g+i,M-i {dx) 
is not necessarily a Levy measure corresponding to a variance gamma process. 


Next we calculate the characteristic function (pj-t of L under P*: 
Proposition 4.5. For any t e [0, T] and u e R, ivith ^ := v — ia, we have 


4‘T-t{0 



X exp|(T- 


M 


-,-(l+h){T-t)C r 


1-i 




G-il 


1 - 




M-1 




M-G-2 


(G + 1)(M-1). 


h(T-t)C 


where p* = Jj[^g{x — e^ + l)v^*{dx). 















4 VARIANCE GAMMA MODELS 


19 


Proof. First of all, we have 

a—Mx POO p—{M—a—iv)x _ „—Mx 


pOQ p — iVlX pO 

/ - dx = 

Jo X Jo 


■ dx 


oo p—{M—ot—iv)x _ p—{M—x)x _j_ p—{M—ix)x _ p—Mx 


= / 

Jo 

pOO pV pOO pM 

= i e-(M-x)x / e-t^dtdx 

Jo Jo Jo JM-dl 

PV poo pM poo 

i / / / e-^^dxdt 

Jo Jo Jm-olJo 

log 


dx 


= I 


= log 


M — a 
M — oi — iv 


f M poo 

iM-a Jo 

M 


M — a 


->■>8(1-1)- 
( 1 ) 


which provides 


p pO pGx poo 

/ (e'?" - l)vc,GMidx) = C / -if—dx + C - 1) 

J]Rn J — OO — X 70 


o-Mx 


■ dx 






= -C(^log (^1+cj+log (^1-M 

In addition, we have 

f xvq G M{dx) =—C f e^^dx + c[ e~^^dx = —C 

7Ro ' J— oo Jo 

Together with Proposition |4]3l we obtain 

[ - 1 - i^x)v^* {dx) = log ("l -I ^ 

JRn V C 


M-G 
GM ■ 




‘°di-M 


- log 1 -H 


If 


G-Il 


hC 


log 1 - 


If 


/iC 




GM 


from which Proposition l4.5l follows. 


M- 1 
M-G-2 
(G + 1)(M-1)' 

□ 


Now, we reformulate iflSl into a linear combination of two Fourier trans¬ 
forms in order to allow use of the FFT. As preparation, we show the following: 


Lemma 4.6. 


[ — l)v{dx) = Glog 

Jrq 


M-i^ G + 
M-l-i^G + l + i^ 


( 2 ) 
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Proof. First of all, we have 

/ -l)v{dx) 

Jrq 

= I - l)v(dx) 

jRn 


= c 


[jo X Jo X J 


/O X 

To calculate Q, we compute 

coo gJC _ ^ 


(3) 


/ 

JO 


e cos bx dx and 


/o X 

for a > 1 and 6 ]R. First, we have 
_ 1 


/' 

JO 


e^-1 


e sin l7x dx 


I 


-e ‘'^coshxdx 


0 X 


/■“ cos fox /■“ 
JO X Ja-\ 


xe ^^dtdx = 




pa pc 

Ja-1 Jo 
2 _|_ ^2 


COS fox ■ e dx dt 


{a — 1)2 + fo2 y ■ 


A similar calculation implies that 

roo gX _ 


/ 

JO 


-e “^sinfoxdx = 


X 


/ 

Jfl-i 


-1 t 2 -|- ^2 


fo . _i fl — 1 

dt = tan -— tan —;—. 


(4) 


(5) 


Noting that M — a > 2 and tan ^ x = j log for x G IR, we have, by iSll and 

©, 


/' 

JO 


- 1 




pOO pX _ 1 nCO pJi _ I 

/ - COSvxdx-\-i / - e~^^~ sinvxdx 

Jo X Jo X 


- 1 


^log 


(M — a)^ + 

(M — a — 1)2 + z;2 
M — a — iv 


= log 


_iM—a _iM —a —1 

z tan-tan - 


M — a — 1 — fiz 


( 6 ) 


Calculating the first term of the right-hand side of (|3]l in the same way as the 
above, we obtain 


coo l^^_(G+.+l+Z.).^^ ^ j / G + X + iv \ _ 

Jo X °\G-|-a-|-l-|-m/ 


JO X 

Substituting ||5]l and ijTll for Q, we arrive at lO. 


(7) 

□ 
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From the above lemma, I 2 is given as follows: 

I Ep. [(Sre^ - K)+- {St - K)+\ Et-]{e’^ - l)v{dx) 

VRo 

1 poo p 

= L {E^’^-l){e’‘-l)v{dx)xp2{0dv 

7T JQ J ]Rq 

= ((AT^wny) 

( 8 ) 


where 


^vciO ■= Clog 


/ M-iC G + fg \ 
'vM-l-z^G + l + z^j 


MO- 


Recall that tp2{0 = ^ result, we need only use the FFT twice for 

computing 12 - 

As the final item of this subsection, we estimate a sufficient length for the 
integration interval of lO for a given allowable error e > 0 in the sense of 0. 
We first provide an upper estimate of (pj-t follows: 


Proposition 4.7. For any iz e IR, 


|^T-t(r’-ia)| < C2|iz| 


where 


X exp I (T — t)ci 


* , M-G-2 

y +{l + h)C - hC- 


GM 


■(G + 1)(M-1)_ 


(9) 


Proof. This can be seen because 


iv + a 


for any a > 0. 

We need to prepare one more lemma: 

Lemma 4.8. 


< 


G“ 


I e'^^{e^-l)v{dx) 

Rn 


< G 


G + a M — a — 1 


□ 


( 10 ) 
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Proof. The same sort of calculations as in ([T) imply 

[ — l)v{dx) 

JlRo 

r 

Jo 


{l/. 

{/o 


1 — 

X 

°° — 1 


< C 

< C 

= C^log(l 
1 


^-{G+a+l+iv)x 


-1 




e“(G+“+i)^dx+ / 

Jo 

log (l + 


< C 


+ 


G “h iT 
1 


1 

X 

1 




M-a-1 


G + a M — a — 1 


□ 


When we calculate lO, N and rj should be taken so that Nt] satisfies (flTTl below 
for a given allowable error e > 0. 

Proposition 4.9. Let e > 0. When a > 0 satisfies 


CC 2 K 


-afi-lgoc 


1 1 

+ TT-T + 


7Te(2C(r-t)+1) [G + a M-a-l 
we have 


log 


MG 


(M-1)(G + 1) 


< fl2C(r-o+i 


1 rOO p 

- ie^^^-iMe^-l)v{dx)^p2{Qdv 

71 J a ^IRq 


< £, 


(11) 


( 12 ) 


where C 2 is defined in dPI. 
Proof. By dlOll . we have 
1 


r / (d?" - l)(e^- - l)vidx)il,2{fi)d 

J a Ro 

-l)v{dx)ip 2 {fi)d' 

J a ^Rn 


+ 


< 1 

71 


f K / {e^-l)v{dx)ip2{0d' 

J a ^Ro 

[ — l)v(dx) + [ (e’^ — l)v{dx) 

JRo .'Ko 

log 


\K 




< - { 
71 


1 ^ 1 
G + a M-a-1 ^ 


MG 


(M-1)(G + 1) 


iJ’iiOldv 

pOO 

/ \uo\ 

J a 


dv I . 

( 13 ) 


Because Proposition l4.7l implies 


\uo\ = 


<pT-t{Qsf. 




< _Lc2|i,|-2C(r-f)sf_ = C2 SMi^|-2G(^-0-2, 
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we 


have, together with ifTSll , 


R.H.S.ofC2) < -CC 2 K 


-a+1 


G + a M — a — 1 

r \^\-2C(T-t)-2dv 
J a 


-CC2R^*+^Sf_ 


G + a M — a — 1 


+ 


+ 


log 


log 


MG 


(M-1)(G + 1) 


MG 

(M-1)(G + 1) 


'2C{T-t) + l' 


□ 


4.2 Numerical results 

We illustrate our numerical results for a variance gamma model. Choosing the 
model parameters as k = 0.15, m = —0.2, and S = 0.45, which meet the second 
condition of Assumption ! 1.11 we compute LRMt for the same numerical exper¬ 
iments as in Subsection 3.2. Note that M > 4 is satisfied. Moreover, we also 
fake the same parameters related to the FFT as in Subsection 3.2. Nt] satisfies 
im for any parameter set. The results are shown in Fig. |2l The computation 
time to obtain Fig.|2lb) was 0.19 s. 

In addition, we implemented the same type of numerical experiments as 
the above based on market data. We used the Nikkei 225 index for March 
2014. We need to set the log price It := log(St/S q), where Sq is the price on 28 
February 2014, which was 14841.07. We estimate the parameters C, G, and M 
in Table [ 1 ] from the mean, variance, and skewness of the log price by using the 
generalized method of moments and the Levenberg-Marquardt method. 


Table 1: Estimated parameters 


C 2.469395026815120 
G 23.743109051760964 
M 24.903251787154687 


Because G — Mw— 1.16, this parameter set satisfies Assumption 11.11 We take 
T = 1 and St- = 14841.07, that is, Lt- = 0. First, fixing the strike price 
K = 14000, we compute LRMt for f = 0,0.05,..., 0.95. Nexf, fixing t to 0.5, the 
values of LRM 0.5 are computed for K = 10000,11000,.. .,20000. Note that Nrj 
satisfies (fTTIl . The results of the computation are illustrated in Fig.|3] 
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(a) Results of the same numerical experiments as Fig.^a) for a variance gamma model 
with parameters k = 0.15, m = —0.2, and S = 0.45. 
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(b) Results for the computation as Fig.^h) for the same variance gamma model as in (a). 
Figure 2: Variance gamma model with parameters k = 0.15, m = —0.2,S = 0.45 
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0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

t 


(a) Values of LRMt with strike price K = 14000 and St- — 14841.07 for t 
0,0.05,..., 0.95. The values of the three parameters C, G, and M are given in Tabled 
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(b) LRMo .5 for K = 10000,11000,... ,20000 with Sq.s = 14841.07. 


Figure 3: Variance gamma model based on the Nikkei 225 index for March 2014 
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